Dispersion relation for anisotropic media 
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The electromagnetic wave propagation in an anisotropic dielectric media with two generic matrices 
e 13 and /j, 13 of permittivity and permeability is studied. These matrices are not required to be 
symmetric, positive definite, and even invertible. In the framework of a metric-free electrodynamics 
approach, compact tensorial dispersion relation is derived. The resulted formula are useful for a 
theoretical study of electromagnetic wave propagation in a classical media and in a modern type of 
media with a generic linear constitutive relation including metamaterials. 
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I. INTRODUCTION 

Development of the modern microscopic technology 
(nano-technology) provides a possibility to manufacture 
materials of a rather non-ordinary electromagnetic pa- 
rameters. This situation recall the theoretical investi- 
gation of wave propagation in a media with a generic 
constitutive law. 

A wide class of media is characterized by a linear con- 
stitutive law and in general can be described by four 3x3 
matrices. Two of these matrices, e lJ and (jf J (permittiv- 
ity and permeability matrices), describe the pure elec- 
tric and the magnetic properties of the matter. Two 
additional matrices describe relativity smaller electric- 
magnetic cross-term effects. 

The ordinary textbook's description of a media with 
two anisotropic matrices e u and // J is based on a diago- 
nalization of one of them [1], [2], [3]. This algebraic pro- 
cedure is always possible for a symmetric matrix. More- 
over, if both matrices are symmetric and one of them is 
positive definite, both of them can be diagonalized. Even 
with this simplification, the corresponded dispersion re- 
lation is given in a rather complicated form. Moreover, 
it is clear that the diagonalization technique is not ap- 
plicable in a general case when both matrices e 1 ^ and 
are not symmetric nor positive definite. 

In the current paper, we study the wave propagation 
in a generic media in the framework of premetric elec- 
trodynamics approach [4], [5], [9]. Our final result is a 
compact form of the dispersion relation. For two generic 
matrices e lJ and fj, lJ , it is given by the expression 

- 2{^k i k j )w 2 + tJhUL ^ ^ = o , (1.1) 



where 



det e det /i 



- imn jpq -1-1 
) nq r^mp ' 



(1.2) 



The organization of the paper is as follows: In the next 
section, the metric-free electrodynamics notations is re- 
called. In Section 3, the covariant metric- free form of 



the dispersion relation is represented. The main results 
are given in Section 4 where several compact forms of 
the generic dispersion relation are derived. In Section 
5, the examples for isotropic, diagonal anisotropic, and 
non-diagonal (magnetized ferrite) media are represented. 



II. ANISOTROPIC MEDIA IN THE 
METRIC-FREE DESCRIPTION 

Let us start with a metric-free four dimensional system 
of Maxwell equations 



0, 



H 



a/3 







(2.1) 



It includes two antisymmetric tensors — the field strength 
J- a /3 and the field excitation Ti a @ . The Greek indices 
change in the range a,/?, ••• = 0,1,2,3, the comma 
denotes the partial derivatives relative to the coordi- 
nates {x°, x 1 , x 2 , x 3 } = {ct,x,y,z\. In sequel, the Ro- 
man indices will be used for the spatial coordinates, 
i,3,"- = 1,2,3. 

The (1 + 3)-decomposition of the field tensors reads 



E{ = T§i , B l - 

D l = H 0i , Hi = 

The electric current is given by 

j° = p, r 



i 



c 



(2.2) 
(2.3) 

(2.4) 



In this notation, the system (2.1) is rewritten in the or- 
dinary three dimensional form of Maxwell equations 



div B = , 
div D = , 



1 <9B 

cur/E+-— = 0, (2.5) 
c at 

curlH-if^j. (2.6) 
c at c 



For a dielectric media, two antisymmetric tensor fields 
are assumed to be linearly related one to another 



1 



(2.7) 
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The constitutive tensor is antisymmetric in two 

pairs of indices, so it has, in general, 36 independent 
components. Such generic constitutive tensor can be rep- 
resented by four 3-dimensional matrices of 9 independent 
components. We will use a representation of the form 

^=(^ S)- (2 - 8) 

In this paper, we restrict to an electromagnetic media 
which describes by two tensors s lJ and 7Ty. Two addi- 
tional tensors j l j and j l j represent the electric-magnetic 
cross-terms, which are relativity small for most types of 
the dielectric materials. We will consider, however, a 
some type of a generalized anisotropic media. In partic- 
ular, we will not require the matrices e y and tt^ to be 
symmetric, positive definite, nor even invertible. Conse- 
quently we will use a constitutive tensor of 18 indepen- 
dent components 

X° i0j = £ ij , X ljkl = -e ijm e kln n mn • (2.9) 

In three-dimensional form, the corresponded constitutive 
relation is given by 

D'^e^Ej, if* = TTyB*. (2.10) 

For a regular matrix 7Ty , an inverse permeability matrix 

(■K- X ) ij = tfi (2.11) 

is defined. With this notation, the constitutive relation 
takes the ordinary form 

D l = e ij Ej , B l = n^Hj . (2.12) 

III. A GENERAL DISPERSION RELATION 

A covariant dispersion relation for a generic constitu- 
tive tensor x a/ljS recently accept a considerable interest 
[5] , [6] . Here we briefly recall the necessary notations and 
the main stages of the derivation as it given in [9]. 

Our aim is to establish the necessary conditions for 
existence of physically non-trivial solutions of the source- 
free system 

e Q ^ 7 % 7 ,5 = , x af3lS F^.s = • (3.1) 

Here the ordinary condition of the geometric optics ap- 
proximation is accepted. In particular, we consider the 
media parameters encoded in x a ^ 1& as varied slowly rel- 
ative to the change of the electromagnetic field. 

The first equation of (3.1) has a standard solution in 
term of the vector potential A a 

Fa/3 = ^ (A*,/3 - -4/3, a) ■ (3-2) 

Consequently, the second equation of (3.1) takes the form 
X^Ajjs^O. (3.3) 



Let us look for a solution of this equation in the form of 
a monochromatic wave ansatz 

A a = a a e iq ^ . (3.4) 

We substitute this ansatz into (3.3) and treat the am- 
plitude of the field a a and the wave covector qp as slow 
functions of a spacetime point. Consequently, we come 
to an algebraic system 

M aS a s = (3.5) 

with a characteristic matrix 

M as = x ^ qpqi . (3 . 6) 

This matrix evidently satisfies the relations 

M aS q a = , M aS q 5 = . (3.7) 

These relations have a clear physical meaning. The first 
equation represents the charge conservation law, while 
the second one means that an ansatz (3.4) with q a ~ a a 
is a solution of (3.5). Certainly this solution is not phys- 
ically meaningful, because it corresponds to a zero value 
of the field T a [i, i.e., it is related to the gauge invariance 
of the field equations. 

Thus we are looking for a solutions of the system (3.5) 
constrained by the relations (3.6). On the matrix lan- 
guage, these relations mean that the columns and the 
rows of the matrix M aS are linearly dependent, i.e., the 
matrix is singular. Consequently, a system always has 
a non-zero solution. However we need more of that, in 
fact, we are looking for an additional linear independent 
solution. Only this one will be of a physical meaning. 

It is an algebraic fact, that a linear system has two 
independent solution only if the adjoint of the character- 
istic matrix equal to zero. So we come to an equation 

(adj M) a0 = . (3.8) 

On a first view, it seems that we require here 16 condi- 
tions for 16 components of the matrix M a/3 . In fact, the 
situation is much different. It can be proved [9] that, for 
a matrix which satisfies the conditions (3.7), the adjoint 
matrix is of the form 

(adj M) a/3 = X(q)q a q /3 . (3.9) 

Consequently, a necessary condition for existence of a 
physically meaningful solution for a wave propagational 
system is expressed by a scalar equation 

\{q)=0. (3.10) 

The function X(q) is a homogeneous 4-th order polyno- 
mial in the wave covector q a . It's explicit forms are given 
in [5], [9]. 
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IV. AN ANISOTROPIC DISPERSION 
RELATION 

For a generalized anisotropic media with a constitutive 
tensor (2.9), we apply an ordinary (l+3)-decomposition 
of the wave covector 

q a = {w,k l ). (4.1) 

The characteristic matrix has now the entries 

M 00 = e^kikj , M° l = -e l3 k 3 w , ikf = -e ij k jW , 

(4.2) 

and 

M« = -e ij w 2 + e mn ei pq ir mp k n k q . (4.3) 
We write the latter equation in a short form 

M ij = e ij w 2 - t 11 . (4.4) 
where a matrix r y is defined as 

T ij = ^mn^p^^ (4 5) 

Due to the relations T y 'fcj = r'-'fej = it is singular. 

In correspondence with (3.9), it is enough to calculate 
only one component of the adjoint matrix. Write 

(adj M) 00 = ^jCn j 2 i 3 £ji j 2 j 3 M iljl M* 2 - 72 M* 3 - 73 . (4.6) 

Substituting (4.4) we derive for (adj M) o = Au> 2 

\ — i-,/, 4 (c. ■ ■ f . . . p i\h p ilh P izh\ _ 
i_,,,2 . . , . . . F ilh F i2h T i333\ i 

In the first term, we recognize the determinant of the 
matrix e. In the second and the third terms, the adjoint 
of the matrices e and r are extracted. Consequently, 
the desired dispersion relation obtains a compact matrix 
form (the central dot denotes the matrix multiplication) 

w 4 (det e) - w 2 tr (t t ■ adj e) + tr (e T ■ adj r) = (4.8) 

This equation can be rewritten in the explicit tensorial 
form. The adjoint of the matrix r is calculated by sub- 
stituting (4.5) into the definition 

(adj t)h = ^e ii2i3 e, j2j3 T^V» 3 . (4.9) 

The result is a scalar function multiplied by kikj 

(adj T )ij = {(adj T,) mn k m k n )kikj . (4.10) 

For a invertible matrices e, fi, we use (4.10) to rewrite the 
expression (4.8) in a form 



where 



— 1 »mn jpq -1 -1 
V — 2 b nqh i mp ■ 



(4.12) 



Note some straightforward facts resulted from this ex- 
pression: 

(1) The dispersion relation (4.11) is symmetric under 
interchange between e and [i. 

(2) A criterion for the absence of zero-frequency modes 
with w = for some (ki,k 2 , ks) ^ (0, 0, 0) takes the form: 
The symmetric parts of the matrices e and \i are definite 
(positive or negative) .[10] 

(3) The necessary and sufficient condition for hyperbol- 
icity (four real w for any real k) is expressed as a system 
of inequalities 



ip zj kikj > positive definite (4-13) 



and 



„ pijhh- n mn h h 

U^kikj) 2 > £ -P^ g , m " > . (4.14) 
yr JJ ~ dete detfj, y 1 

(4) In the non-birefringence case [11] with a unique 
optical metric, the dispersion relation (4.11) takes the 
form 



w 



2 hkj = 0. 



(4.15) 



Thus the optical metric is of Minkowski signature if and 
only if the matrix ip l i is positive definite. 



V. EXAMPLES 

A. Isotropic case 

It easy to check that in the isotropic case with = 
sfiij — tio~ij (4-11) yields the ordinary dispersion re- 
lation (efi)w 2 - k 2 = . 



B. Diagonal case 

Let us consider a more involved example of two diago- 
nal matrices 

e = diag(e 1 ,e 2 ,e 3 ) , fj, = diag(m, fj, 2 , Ate) • (5-1) 
The last term of (4.11) takes the form 

£\k\ + e 2 k\ + e 3 k 2 \i x k\ + [i 2 k\ + fi 3 k 2 



Si£ 2 £3 



A*iAteAte 



(5.2) 



The coefficient of the second term of (4.11) also easily 
calculated 



V^Ate £3^2/ V £ iAte £3Ml/ \ £ 1/- 1 2 £2t 



£2^1/ 
(5.3) 
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Consequently, the dispersion relation in this case takes 
the form 

(ei£2£3MiA t 2M3)w 4 + (sifnkl(e 2 lJ-3 + e 3 fi 2 ) + 

£2^1 ( £ iM3 +£3^1) +£3M3fcI( £ iM2 +£2A*i)) w2 + 
(eifcf +£ 2 fc| +£3^1)^1^ + M2fc| + Ms*!) =0 (5.4) 

In a special case m = l, this formula coincides with the 
one given in [2], [3]. 

C. Magnetized ferrite 

Magnetized ferrite materials are described by an 
isotropic dielectric constant 



e = diag{e 1 e, e) . 



(5.5) 



Under influence of the magnetic field the initially 
isotropic magnetic matrix obtains an anisotropic mod- 
ification. For a magnetic field directed as the z-axis [12], 



jU iq 
— iq fi 
^0 



(5.6) 



The last term of (4.11) takes now the form 

P_ fl(k 2 + fc|) + flgk 2 

e 2 fi (^ 2 - <7 2 ) 



(5.7) 



The coefficient of the second term of (4.11) is easily cal- 
culated. The inverse matrix 



r MoM I 



H iq 

-iq fx I where r = 



fJ-o 



(5.8) 



Consequently the second term of (4.11) takes the form 



( S mp S n q _ S m H np )—^ p k n k q 



W 



K k l + k l) +M0&3 



The resulting dispersion relation is 



(5.9) 



r^ 2 s 2 w 4 - w 2 e^ [2k 2 + (r - l)(fc 2 - k 2 )] + 

'Mo 
H 



k 2 



k 2 +[^-l)k 2 



0. 



(5.10) 



In absent of an exterior magnetic field, q — and \i — /i . 
Consequently, r = 1 and the isotropic dispersion relation 
reinstated. 

VI. CONCLUSION 



A short dispersion relation is derived for a general- 
ized anisotropic media. The corresponded matrices are 
not required to be symmetric, positive definite, and even 
invertible. This compact form is much simpler for cal- 
culation and for theoretical analysis than the ones rep- 
resented in the electromagnetic literature [13], [14]. The 
algebraic consideration presented here can be also use- 
ful for a compact algebraic representation of a general 
bianisotropic tensorial dispersion relation represented re- 
cently in [15]. 



Acknowledgments 



I thank F.-W. Hehl, V. Perlick, Y. Obukhov and Y. 
Friedman for most useful discussions. 



[1] A. Sommerfeld , Electrodynamics (New York: Academic 
Press, 1964). 

[2] L. D. Landau and E. M. Lifshitz, Electrodynamics of 
Continuous Media (Pergamon, Oxford, 1984). 

[3] Jackson J. D. Classical electrodynamics(3ed., Wiley, 1999) 

[4] E.J. Post, Formal Structure of Electromagnetics - Gen- 
eral Covariance and Electromagnetics (North Holland: 
Amsterdam, 1962). 

[5] F.W. Hehl and Yu.N. Obukhov, Foundations of Classical 
Electrodynamics: Charge, Flux, and Metric (Birkhauser: 
Boston, MA, 2003). 

[6] Y. N. Obukhov, T. Fukui and G. F. Rubilar, Phys. Rev. 
D 62, 044050 (2000) [arXiv:gr-qc/0005018]. 

[7] F.W. Hehl and Y.N. Obukhov, Phys. Lett. A334 (2005) 
249-259; arXiv.org/physics/0411038. 

[8] Y. Itin, J. Phys. A 40, F737 (2007) [arXiv:0705.0756 



[hep-th]]. 

[9] Y. Itin, "On light propagation in premetric electrody- 
namics. Covariant dispersion relation," arXiv:0903.5520 
[hep-th]. 

[10] V. Perlick, Privite communication. 

[11] Y. Itin, Phys. Rev. D 72, 087502 (2005) [arXiv:hep- 
th/0508144]. 

[12] A. Ishimaru, Electromagnetic Wave Propagation, Radia- 
tion, and Scattering ( Prentice Hall, New Jersey, 1991). 

[13] Damaskos, N. J.; Maffett, A. L.; Uslenghi, P. L. E., IEEE 
Trans. Antennas Propagat., AP-30 (1982), 991-993. 

[14] Graglia, R. D.; Uslenghi, P. L. E.; Zich, R. E., IEEE 
Trans. Antennas Propagat. 39, (1991), 83-90. 

[15] Y. N. Obukhov and F. W. Hehl, Phys. Rev. D 70, 125015 
(2004) [arXiv:physics/0409155]. 



